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Abstract 



Using plaquette and Symanzik improved gauge action and stout link clover 
fermions we determine the improvement coefficient csw in one-loop lattice per- 
^ ' turbation theory from the off-shell quark-quark-gluon three-point function. In ad- 

dition, we compute the coefficients needed for the most general form of quark field 
improvement and present the one- loop result for the critical hopping parameter Kc- 
We discuss mean field improvement for csw and Kc and the choice of the mean field 
coupling for the actions we have considered. 



1 Introduction 



Simulations of Wilson-type fermions at realistic quark masses require an improved action 
with good chiral properties and scaling behavior. A systematic improvement scheme that 
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removes discretization errors order by order in the lattice spacing a has been proposed 
by Symanzik [1] and developed for on-shell quantities in [2l[3]. 0{a) improvement of the 
Wilson fermion action is achieved by complementing it with the so-called clover term [3], 
provided the associated clover coefficient is tuned properly. 

Wilson-type fermions break all chiral symmetries. This introduces an additive negative 
mass renormalization term in the action, which gives rise to singularities in the quark 
propagator at small quark masses and makes the approach to the chiral regime difficult. 
A chiral improvement of the action is expected to reduce the additive mass renormalization 
and the spread of negative eigenvalues. Surprisingly, this is not accomplished by the clover 
action. 

While the magnitude of the additive mass term decreases with increasing clover term, 
the problem of negative eigenvalues is more severe for the clover than for the standard 
Wilson action. It is well known that via a combination of link fattening and tuning of the 
clover coefficient, it is possible to reduce both the negative mass term and the spread of 
negative eigenvalues [H El E] . 

The focus of this investigation is to determine the clover coefficient and the additive 
mass renormalization for plaquette and Symanzik improved gauge action and stout link 
clover fermions in one-loop lattice perturbation theory. 

The Symanzik improved gauge action reads [1] 



< Co ^ - Re Tr (1 - f/piaqucttc) + Cl ^ - Re Tr (1 - f/Rcctangle) \ (1) 
L Plaquette Rectangle ) 



C-Sym _ 6 

g , 

' Plaquette Rectangle 



with Co -|- 8ci = 1 and 
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co=3, c, = --. (2) 

This reduces to the standard plaquette action Sq ^"^ for ci = 0. 
Clover fermions have the action for each quark flavor [3] 

Sf = "^^ X^ I - ^ ^(x)f/^(x) (1 -7/,)^(x-ha/i) 

X 

+ il>{x)Ul{x - afx) {I + '^ ^) il){x - afi) (3) 

+ - (4 -h aniQ + am) ^(x)V^(a;) - csw 9 j ip{x) a^uF^u{x) i){x] 



where 

o-^o = 7, 4 , (4) 

ZKc 

Kc being the critical hopping parameter, is the additive mass renormahzation term, and 
Ffj,u{x) is the field strength tensor in clover form with a^^ = {i/2) (7^7,^ — ■Julu)- We 
consider a version of clover fermions in which we do not smear links in the clover term, 
but the link variables [/^ in the next neighbor terms have been replaced by (uniterated) 
stout hnks [7] 

^^(x) = e*«-(^)f/^(x) (5) 

with 



Q,{x) = 1^ 



V,{x)Ul{x) - f/,(x)V;t(x) - ^Tr {V,ix)Ulix) - U,{x)Vi{x)) 



(6) 



V^(x) denotes the sum over all staples associated with the link and cj is a tunable weight 
factor. Stout smearing is preferred because ([5]) is expandable as a power series in g"^, 
so we can use perturbation theory. Many other forms of smearing do not have this nice 
property. Because both the unit matrix and the 7^ terms are smeared, each link is still a 
projection operator in the Dirac spin index. 

The reason for not smearing the clover term is that we want to keep the physical 
extent in lattice units of the fermion matrix small which is relevant for non-perturbative 
calculations. In that respect we refer to these fermions as SLiNC fermions, from the 
phrase Stout LinkNon-perturbative Clover. 

The improvement coefficient Csw ^^ well as the additive mass renormahzation aniQ 
are associated with the chiral limit. So we will carry out the calculations for massless 
quarks, which simplifies things, though it means that we cannot present values for the 
mass dependent corrections. 

For complete 0{a) improvement of the action there are five terms which would have 
to be added to the 0{a) effective action, they are listed, for example, in [8]. Fortunately, 
in the massless case only two remain, 

Ci = ^a^.F^.^, (7) 

O2 = iJDD^p. (8) 

The first is the clover term, the second is the Wilson mass term. We have both in our 
action, there is no need to add any other terms to the action. 

In perturbation theory 

csw = l + 9'4w + Oig^). (9) 



The one-loop coefficient Cg{y has been computed for the plaquette action using twisted 
antiperiodic boundary conditions [9] and Schrodinger functional methods [10]. Moreover, 
using conventional perturbation theory, Aoki and Kuramashi [TT] have computed c^l^ for 
certain improved gauge actions. All calculations were performed for non-smeared links 
and hmited to on-shell quantities. 

We extend previous calculations of c^^y to include stout links. This is done by comput- 
ing the one-loop correction to the off-shell quark-quark-gluon three-point function. The 
improvement of the action is not sufficient to remove discretization errors from Green 
functions. To achieve this, one must also improve the quark fields. The most general 
form consistent with BRST symmetry is [12)H 

V^^(a;) = [l + acD IP +aig cngiA{x) ] tpix) . (10) 



From now we denote improved quark fields and improved Green functions by an index -k. 
These are made free of 0{a) effects by fixing the relevant improvement coefficients. 

There is no a priori reason that the gauge variant contribution cngi4-{x) vanishes. 
The perturbative expansion of cngi has to start with the one-loop contribution [12]. As 
a byproduct of our calculation we determine that coefficient cj^qj 

CNGI = g' c^NGi + Oig') (11) 

and find that it is indeed nonvanishing. 



2 Off- shell improvement 

It is known [H] that the one-loop contribution of the Sheikoleslami-Wohlert coefficient in 
conventional perturbation theory can be determined using the quark-quark-gluon vertex 
f^^{pi,P2,csw) sandwiched between on-shell quark states, pi (^2) denotes the incoming 
(outgoing) quark momentum. In general that vertex is an amputated three-point Green 
function. 

Let us look at the 0{a) expansion of tree-level Ajj, {pi,p2,csw) which is derived from 
action ([3]) 

Aa\pi,P2,csw) = -ig'^^,- g\al{pi +P2)^, + cswig^aa^aipi -P2)a + 0{a^) . (12) 



^In [12] the authors use and instead of Ip and 4 - both choices are equivalent. Our choice is 
motivated by the discussion of off-sheh improvement in the next section. 



For simplicity we omit in all three-point Green functions the common overall color matrix 
T". That tree-level expression between on-shell quark states is free of order 0{a) if the 
expansion of csw starts with one, as indicated in Qi 

u{p2) ^%{pi,P2)u{pi) =u{p2){-ig-fi,)u{pi). (13) 

Therefore, at least a one-loop calculation of the A^(pi,p2, c^^-) is needed as necessary 
condition to determine Csiy- 

The off-shell improvement condition states that the non-amputated improved quark- 
quark-gluon Green function G^,f^{pi,p2,q) has to be free of 0{a) terms in one-loop accu- 
racy. In position space that non-amputated improved quark-quark-gluon Green functions 
is defined via expectation values of improved quark fields and gauge fields as 

G,^{x, y, z) = (^,(x) ^M M^)) ■ (14) 

Since the gluon propagator is (9(a)-improved already, we do not need to improve gauge 
fields. Using relation (TTOl) we can express the function G^^ by the unimproved quark fields 



G.^{x,y,z) = G^{x,y,z) + acD((^l/)l/) ' + IP 'i/))a, 

+ iag CNGi ( (4^"' + ^~ V) Am) . (15) 

where Gn{x, y, z) is the unimproved Green function. 
Taking into account 

' iiU)'^ + ^" V) A^) = 2 a CD 5{x - y) {A^{z)) (16) 



and setting {A^{z)) = (unless there is an unexpected symmetry breaking), we obtain 
the following relation between the improved and unimproved Green function 

G^^ix, y, z) = G^{x, y,z) + iag cngi ( (4^"' + 0" V) A^) • (17) 

From fITTj) it is obvious that tuning only csw to its optimal value in G^{x,y,z), there 
would be an 0{a) contribution left in the improved Green function. The requirement 
that G^^{x,y, z) should be free of 0{a) terms leads to an additional condition which 
determines the constant cngi- It has not been calculated before. 



Taking into account the expansion (ITT]) of cngi we get in momentum space (JF[- 
denotes the Fourier transform) 



i a g Cngi ^ 



(4^-^ + ^-^)^^^] =^«^'c!iG/(7.-^ + T^7.) KX\q). (18) 



or its amputated version 



i a g Cngi ^ 



(4^"' + ^^ V) A^)'JJ =-a/cii^, (^2 7/^ + 7mI*i) • (19) 



The relation between non-amputated and amputated unimproved and improved three- 
point Green functions are defined by 

Gf,{pi,P2,q) = S{p2)K{pi,P2,q,c^sw)S{Pi)Ky^{q), (20) 

G^^{pi,P2,q) = S4p2)A^,yipuP2,q)S^{pi)K^^{q), (21) 

Ky^{q) denotes the full gluon propagator which is 0(a)-improved already, S{p) and Si,{;p) 
the corresponding quark propagators. 

With the definition of the quark self energy 

S(p) = iso + z^Si(p) + '^Mp) (22) 

the unimproved and improved inverse quark propagators are given by 

S^\p) = ^t^lip) + ^^2{p) = ^f^lip)l^l-la^f^y (23) 

S:\p) = ifT.^{p). (24) 

Using the Fourier transformed fll7j) with fll9p and amputating the Green function (1201) . 
taking into account the inverse quark propagators (!23l) . we get the off-shell improvement 
condition in momentum space 

AM(Pl'P2,g,4w) = A^M(Pl'P2,g) + a^^cS^G/(|*2 7M + 7M^l) 

- 2 * ^2 ^ r \ ^*m(pi , P2, g) - 2 A*m(pi , P2, g) « ^1 slvT ■ 

This expression should hold to order 0{g^) by determining both Cj^qj and Cg{y correctly. 
It is clear from (125!) that the improvement term oc Cnqj does not contribute if both quarks 
are on-shell. 



3 The one-loop lattice quark-quark-gluon vertex 



The diagrams contributing to the amputated one-loop three-point function are shown in 
Fig. [H The calculation is performed with a mixture of symbolic and numerical tech- 



q=P2-Pi 







(d) 





Figure 1: One- loop diagrams contributing to the amputated quark-quark-gluon vertex. 



niques. For the symbolic computation we use a Mathematica package that we developed 
for one-loop calculations in lattice perturbation theory (for a more detailed description 
see [S]). It is based on an algorithm of Kawai et al. [I5]. The symbolic treatment has 
several advantages: one can extract the infrared singularities exactly and the results are 
given as functions of lattice integrals which can be determined with high precision. The 
disadvantage consists in very large expressions especially for the problem under consider- 
ation. In the symbolic method the divergences are isolated by differentiation with respect 
to external momenta. 



Looking at the general analytic form of the gluon propagator for improved gauge 



actions [16] one easily recognizes that a huge analytic expression would arise. As discussed 
in [16] we split the full gluon propagator D^l^{k,^) 

Djr(^,0 = Dl'r{k,0 + AD,^{k) , (26) 

where ^ is the covariant gauge parameter (,^ = corresponds to the Feynman gauge). 
The diagrams with D^l^^{k,^) only contain the logarithmic parts and are treated with 
our Mathematica package. The diagrams with at least one AD^y{k) are infrared finite 
and can be determined safely with numerical methods. The decomposition fl26|l means 
that we always need to calculate the plaquette action result, as part of the calculation for 
the improved gauge action. Therefore, we will give the results for both plaquette gauge 
action and Symanzik improved gauge action using the corresponding gluon propagators 
D^^'^ and -D^^™, respectively. 

Because the numerical part determines the accuracy of the total result we discuss it 
in more detail. There are several possibilities to combine the various contributions of the 
one-loop diagrams as given in Fig. [H In view of a later analysis we have decided to group 
all coefficients in front of the independent color factors Cp and Nc. and the powers of the 
stout parameter uj 

^num. ^ ^^ ^^^ j^ Ciu + C2UJ^ + C^ u"") + A^, (C4 + C^u + C^^uj^ + C^ u^) , (27) 

where the Ci have to be computed numerically. In order to obtain Ci we first add all 
contributions of the diagrams shown in ([1]) and integrate afterwards. We have used a 
Gauss-Legendre integration algorithm in four dimensions (for a description of the method 
see [H]) and have chosen a sequence of small external momenta (pi,P2) to perform an 
extrapolation to vanishing momenta. 

Let us illustrate this by an example: the calculation of the coefficient C4. We know 
the general structure of the one-loop amputated three-point function as (we set a = 1) 

Mf,{pi,P2) = 'yf,A{pl,p2) + lpl,f,B{pl,p2) + lp2,^,C{pl,P2) 

+ (Tf,aPl,aD{pi,p2) + (Tf,a P2,a E{pi, P2) . (28) 

From this we can extract the coefficients by the following projections 

Ti-ff^Mf, = 4:A{pi,p2), /i fixed, 

TrM^ = 4pi,^5(pi,P2)+4p2,mC(Pi,P2), (29) 

^Tra.^M^ = 12p^^, D{pi,p2) + I2p2,, E{p^,p2) . 



Relations (!30|) show that one has to compute the three-point function for all four values 
of fi. Further they suggest choosing the external momenta orthogonal to each other: 
Pi ■ P2 = 0. A simple choice is pi,^ = (0, 0, 0,^1,4) and p2,/i = (0, 0,p2,3, 0). 

We discuss the determination of B{pi,p2) and C{pi,p2) in more detail. For small 
momenta they can be described by the ansatz 

B{pi,p2) = Bo + BipI + B2pI, 

C{pi,P2) = Co + C,pI + C2pI. (30) 

The choice of the momenta is arbitrary except for two points. First, they should be 
sufficiently small in order to justify ansatz (l30l) . Second, they should not be integer 
multiples of each other in order to avoid accidental symmetric results. The symmetry 
of the problem demands the relation Bq = Cq which must result from the numerical 
integration also. Performing the integration at fixed pi and p2 we obtain complex 4x4 
matrices for 71^3(^1,^2) and M4^{pi,p2) from which the quantities B{pi,p2) and C{pi,p2) 
are extracted via fl5U]) . 



A nonlinear regression fit with ansatz ( 13011 gives 

Bo = 0.00553791 with fit error 5Bo = 7 x 10~* , 

Co = 0.00553789 with fit error 5Co = 6 x 10"^ . (31) 

It shows that the symmetry is fulfilled up to an error of (9(10"'') which sets one scale for the 
overall error of our numerical calculations. In Fig. [2] we show the almost linear dependence 
of B{pi,p2) and C{pi,p2) on pf. (In the integration we haven chosen pi^^ = 0.87^2,^* so 
that we can restrict the plot to one variable.) 

Another source of errors is the numerical Gauss-Legendre integration routine itself. 
We have chosen a sequence of n^ = 14^, 18^, 22^, 26^ and 30^ nodes in the four- dimensional 
hypercube and have performed an extrapolation to infinite nodes with an 1/n^ fit ansatz. 
Both procedures, Gauss-Legendre integration and the fit p -^ 0, give a combined final 
error of 10~^. 

The third error source are the errors of the lattice integrals of our Mathematica calcu- 
lation for the terms containing the plaquette propagator D^j^^'^ only. These integrals have 
been calculated up to a precision of (9(10~^°). Therefore, their errors can be neglected in 
comparison with the others. 

Summarizing we find that the error of our numerical procedure is of (9(10~^). Addi- 
tionally, we have checked our results by an independent code which completely numerically 
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Figure 2: B(pi,p2) (circles) and C{pi,p2) (squares) as function oi pf together with their 
corresponding linear fits in p^. 

computes the one-loop contributions for each diagram including the infrared logarithms. 
Both methods agree within errors. 

The Feynman rules for non-smeared Symanzik gauge action have been summarized 
in [H] . For the stout smeared gauge links in the clover action the rules restricted to equal 
initial and final quark momenta are given in [6]. As mentioned in the introduction we 
perform a one-level smearing of the Wilson part in the clover action. The corresponding 
Feynman rules needed for the one-loop quark-quark-gluon vertex are much more compli- 
cated than those in [6]. The qqgg- vertex needed in diagrams (c) and (d) of Fig. [1] receives 
an additional antisymmetric piece. The qqggg-vertex in diagram (e) does not even exist 
in the forward case. The Feynman rules are given in Appendix A. The diagrams which 
are needed for the calculation of the quark propagator are shown in Fig. [31 We have 
performed our calculation in general covariant gauge. 
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Figure 3: One- loop diagrams contributing to the quark self energy. 

4 Results for the improvement coefficients and criti- 
cal hopping parameter 

The anticipated general structure for the amputated three-point function in one-loop is 

CL Q f \ Qj Q 

+ ^'^' 2 16^ 1^2 7,. + 7m ^ J + Ciat ^ 2 16^^ ^^" ^° • ^^^^ 

A^'^(pi,P25 q) is the universal part of the three-point function, independent of the chosen 
gauge action, computed in the MS'-scheme 

K^^{PuP2,q) = -i9'yt^-9-'i-{Pl,^^+P2,t^)-Cswi9-(r^^aqa 



-'- 9 \MS/^ „ ^\ , ^ 9 .A/5 



We have calculated the complete expressions for Af'J{pi,p2,q) and A2^J{pi,p2,q). 

The C(a) contribution, A^^{pi,p2,q), simplifies if we set csw = 1 + Oi^g"^) as in ^. 
After some algebra we find 

Ajf(pi,P2,g) = ^(|*2A5f(pi,P2,g) + A5f(pi,p2,g)^i) 

- Cf {p^ 7m (1 - 0(1 - log(p^//i')) (34) 

+7M^i(l-0(l-log(P?//i'))^ 
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where /i^ is the MS mass scale (not to be confused with the index /i). Therefore, we 
only need Af^(pi,p25?) to present the one-loop result (15^ . ^i^^{pi,P2,(l) is given in 
Appendix B. 

If we insert (132|) and (133|) with (134|) into the off-shell improvement relation (!25|) we get 
the following conditions that all terms of order 0{ag^) have to vanish 

(c!iL/ - ^ (^lat - i?iat - S21)) (^2 7/. + 7m^i) = , (36) 

with S21 defined from fl23l) as 

III - l.^((l-0(l-.og(aV)) + E.,) 

and 

S21 = Cf (-(1 - log(aV) + S2i,o) . (38) 

The constant ^21,0 depends on the chosen lattice action. 

It should be noted that equations (1551) and (15B]) are obtained by using the general struc- 
ture (151|) only - we do not need to insert the complete calculated result for A^(pi, ^2, o)- 
In order to get momentum independent and gauge invariant improvement coefficients we 
see from (!35|) that Ciat itself has to be constant and gauge invariant. From (l36l) and (!38|) we 
further conclude that the log(a^yU^)-terms from Aiat and -Biat have to cancel those from S21. 
The same is true for the corresponding gauge terms. The terms oc (1 — ^)(1 — \og{pf / ji^)) 
{i = 1, 2) coming from (1371) are canceled by the corresponding terms in 



Therefore, the relation between Af^^{pi,p2,q) and A^^{pi,p2,q) as given in (134|) is a 
nontrivial result. Once more, it should be emphasized that this relation only holds if we 
use csw = 1 at leading order in g'^. If ( Bll) were not true, we would not be able to improve 
the Green functions by adding the simple 0{a) terms we have considered. 

For completeness we also give the corresponding one-loop values for the quark field 
improvement coefficient cd as defined in (ITOl) . They can be derived from the 0{a) im- 
provement of the quark propagator. The one-loop improvement coefficient c}^ is related 
to the quark self energy by 

^ f, . 9^^F ^c^ ^ ^^ , a^^_4n _ Wi , 2 (1) 



cz. = -^ (1 + ^ (2Si - S2) j + 0{g') = -- (1 +^^c^^ ) + 0{g') . (39) 
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c^, has been calculated for ordinary clover fermions and plaquette gauge action in [T3] . 

Now we present our numerical results for general covariant gauge C, and as function of 
the stout parameter u. For the plaquette action with stout smearing the quantities Aiat, 
i?iat and Ciat are obtained as 

Af^^'^q = Cp (9.206269 + 3.792010^- 196.44601 CU + 739.683641 cj2 



+ (l-01og(aV) 
+ Nc (-4.301720 + 0.693147^+ (1 - ^/4) log(aV)) , 
B^^P = Cf (9.357942 + 5.727769^ - 208.583208 cu + 711.565256 cu^ 

+2(l-01og(aV)) (40) 

+ iVc (-4.752081 + 0.693147^ + 3.683890 cj + (l-^/4)log(aV)) , 

CjJ'^q = Cf (26.471857+ 170.412296 cj- 582.177099^2) 
+ N^ (2.372649 + 1.518742 cu- 44.971612 cu^j _ 

For the stout smeared Symanzik action we get 

J^y^ = Cf (5.973656 + 3.792010^ - 147.890719 tu + 541.380348 cu^ 

+ (l-01og(aV)) 

+ Nc (-3.08478 + 0.693159 ^ - 0.384236 u + {1 - ^/4) log(aV^)) , 
Sfj^"^ = Cf (6.007320 + 5.727769 ^ - 163.833410 cu + 542.892478 u^ 

+ 2(l-01og(aV)) (41) 

+ iVe (-13.841082 + 0.693179 ^ + 3.039641 u + {1 - ^/4) log(aV^)) , 
Cf^y^^ = CiT (18.347163 +130. 772885 cj- 387.690744 a;2j 

+ Nc (2.175560 + 2.511657CJ- 50.832203 w^) . 

As shown in fl25l) (or equivalently fl36l) ) we need the self energy parts Si(p) and ^2{p) as 
defined in f l23|) to solve the off-shell improvement condition. They have the general form 

Ei(p) = 1 - ^ [(1 - log(ay ) + Si,o] , 

S2(p) = 1 - ^ [2 (1 - log(aV) + S2,o] . (42) 
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For the plaquette and Symanzik actions we obtain 

S^o^^ = 8.206268 - 196.446005 u + 739.683641 uj^ + 4.792010 ^ , 

tF^^P = 7.357942- 208.583208 CU + 711. 565260 cj2 + 7.727769^, 

S^yj" = 4.973689- 147.890720 CU + 541.380518 cj2 + 4.792010^, (43) 

S^yj" = 4.007613- 163.833419 CU + 542.892535 cj2 + 7.727769^. 



..UJ^ 



This resuhs in the following expressions for E21 as defined in fl38|) 

S^j^^ = Cf( -0.151673 -1.935759^ + 12.137203 CJ + 28.118384 ( 

-(l-Olog(aV)), 
eII"^ = C^ I'- 0.033924 -1.935759^+ 15.942699 cj-1.512017cu2 (44) 

-(l-Olog(aV) 



Inserting the corresponding numbers into (!35|) . (136|) and (139|) . we obtain the one- loop 
contributions of the clover improvement coefficient 

cJI^P^^'^ = C^ (0.167635 + 1.079148 cu- 3.697285 cu^) 

+ iVe (0.015025 + 0.009617CJ- 0.284786 cu^) , (45) 

cJJ^Sy"^ = C^ (0.116185 + 0.828129 u;- 2.455080 cu^) 

+ N^ (0.013777 + 0.015905 cu - 0.321899 u^) , (46) 

the off-shell quark field improvement coefficient 

cj^^Pi^'i = iV, (0.001426 -0.01 1664 cj) , (47) 

cS^J;^^"' = iV, (0.002395 -0.010841 cj) , (48) 

and the on-shell quark field improvement coefficient 



c 



c 



(l),Plaq 
D 

(l),Sym 
D 



Cf (0.057339 + 0.011755^ - 1.167149 cu + 4.862163 cj^) , (49) 

Cf (0.037614 + 0.011755^- 0.835571 cj + 3.418757cj2) , (50) 



for the plaquette and Symanzik action, respectively. For u; = both the plaquette result 
f H5|) and the Symanzik result (H6l) agree, within the accuracy of our calculations, with the 
numbers quoted in [9], [10] and [H] . 
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From Ward identity considerations it is known that the coefficient cngi has to be 
proportional to N^ only. Additionally, cjqci and csw should be gauge invariant. Both 
conditions are fulfilled within the errors which have been discussed in the previous section. 
It should be noted that ( H71) and ( HHl) are the first one-loop results for the quark field 
improvement coefficient cngi- The gauge dependent improvement coefficient cd depends 
only on the color factor Cp because it is determined by 0{a) improvement of the quark 
propagator. 

The additive mass renormalization is given by 

This leads to the critical hopping parameter Kc, at which chiral symmetry is approximately 
restored, 

" 8 V 167r2 4 y ^ ^ 

Using the plaquette or Symanzik gauge actions, we obtain 

E^^^"^ = -31.986442 + 566.581765CU- 2235.4070870^2^ (53) 

j.Sym ^ -23.832351 + 418.212508CJ- 1685.597405 cj^ (54) 

This leads to the perturbative expression for Hc 

fi;Piaq = - [1+/Cp (0.050639 -0.896980 CU + 3.697285 cj2)] , (55) 

K^"^ = - [1+/Cp (0.037730 -0.662090 u; + 2.668543 u;2)] . (56) 

For both actions arriQ can be tuned to zero for admissible values of u. Using the smaller 
possible value we find u = 0.089396 for the plaquette action and u = 0.088689 for the 
Symanzik gauge action. 



5 Mean field improvement 

It is well known that one-loop perturbation theory in the bare coupling constant g^ leads 
to a poor approximation. The coefficient of g^ is large in most quantities, and the series 
converges poorly. One traditional way to reduce this problem is by mean field improve- 
ment, which consists of two ideas. 
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The first is that we calculate each quantity in a simple mean field approximation, and 
then re-express the perturbative result as the mean field result multiplied by a perturbative 
correction factor. If the mean field approximation is good, the correction factor will be 
close to 1, and we have resolved the problem of the large one-loop coefficient. As a 
good internal test of this part, we can simply look to see how large the coefficient in 
this correction factor is (the "tadpole improved coefficient"), compared with the initial 
unimproved coefficient. 

The second part of the mean field approximation is that we change our expansion 
parameter from the bare coupling g^ to some "boosted" coupling constant, g\[pi which 
we hope represents physics at a more relevant scale, and leads to a more rapidly con- 
vergent series. A well-chosen boosted coupling would reduce the two-loop coefficient. 
Unfortunately we usually cannot test this part of the improvement procedure, because 
the two-loop coefficient is unknown. Fortunately, if the mean field approximation is good, 
the exact choice of boosted coupling constant will not be too crucial, because the lowest 
order improved coefficient will be a small number. 



5.1 Mean field approximation for smeared fermions 

In the mean field approximation we typically assume that the gauge fields on each link 
are independently fiuctuating variables, and that we can simply represent the links by an 
average value uq. Typical choices for uq would be to choose Uq to be the average plaquette 
value, or to choose uq to be the average link value in the Landau gauge. 

A natural question is how we should extend the mean field approximation if we employ 
smearing. One possibility is to express everything in terms of two quantities, Mq, a mean 
value for the unsmeared link, and M5, a mean value for smeared linkqj- We will discuss the 
relation between these two quantities later, first we want to make a general point about 
mean field approximations and smearing. 

The reason we smear our gauge links is to suppress very short range fiuctuations in 
the gauge field, which is justified by the argument that these short range fiuctuations 
are very lattice-dependent, rather than physical. However, put another way, suppressing 
short range fluctuations means that we are correlating nearby gauge links. So there is a 
certain tension between smearing and the mean field notion that each link is fiuctuating 
independently. We will take the attitude that it does still make sense to use the mean 
field approximation if smearing is mild - but we should treat the results with some degree 



^PR would like to thank Colin Morningstar for conversations on this point. 
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of caution if extreme smearing is used. 

Applying this double-M mean field approximation to the SLiNC fermion matrix we 
find the following results for the principal fermion quantities, 

Ei(p)RiM5, S2(p)^M5, Z^^Us, Hc^- , Csw ^ ^ (57) 

8us Mo 

(we define Z^ by the relation S'''™ = Z^S^^^). For reasonable smearing we expect the 
smeared link us to be closer to 1 than the bare link mq, so most quantities will lie closer to 
their tree-level values with smearing. However, the clover coefficient csw is an exception; 
it will be further from 1 with smearing than without, because we construct our clover 
term from unsmeared links. 

As a result, we obtain the mean field expressions for n^ and csw by performing the 
following replacements 

-| pert/ 2 \ 

«c(/) - Krigl,, us) = I "^^MmA ^^^gl^) (58) 

o Us 

and 

MF/2 „, „, \ "•s' "o yyMF 



us ur'jgj ^^ 

^1 uT\gl,) 



csw{.g ) -^ Cswi9MF^ Us, uo) = -J ' csw{9mf) ■ (59) 



Here us and Uq are the measured smeared and unsmeared links at the given coupling and 
Ug^"^ and Uq'^^^ denote the corresponding expressions in lattice perturbation theory. 



5.2 The smeared plaquette in perturbation theory 

We will use Ug'^^ derived from the smeared perturbative plaquette Ps 



uT ^ pT- (60) 



To one-loop order we have 

vK?'' = 1 - 

167r 

9 4./ Ci Hj 



uT = 1 - ^ ^5 , (61) 



witf§ 



A;5 = 87rV / j—-D^p{k)\Vc,i{k,uj)Vpi{k,uj)si{k) + V^2{k,uj)Vp2{k,uj)si{k) 

- iy^^{k,uj)Vp2{k,uj) + Vp^{k,u)V^2{k,u)) si{k)s2{k)\ (62) 



■^We have written this integral for the case of a plaquette in the 1-2 plane, any orientation gives the 
same result. 
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where Dapik) the gluon propagator for the action in question. The smearing function 
Va^ik, uj) is defined in (lA.SP in Appendix A, s^{k) and s^ik) used below are given in (]A.2p . 
Using symmetry and the definition of V , the expression simphfies to 

ks = levrV j 0- [Dr,{k)s2{k)s2{k) - D,^{k)s,{k)s,{k)] (l - Au s\k)f . (63) 

We can see from this form that mild smearing has the effect of suppressing the contribution 
from large k. Setting cj = in ks-, we recover the unsmeared link in perturbation theory 

ur' = l-^^ks{u = Q). (64) 

For the plaquette action propagator we can calculate the integral exactly. The result 
is 

kl'^'' = Ti^{l-lQuj + 72uj^) . (65) 

Let us see how well this improves the expressions for n^ and csw- Using the result fl55l) 
we find 

^piaq,MF ^ _}_ U^gl^j^Cp (-0.011861 + 0.103020 cj- 0.802715 cj2)l (66) 

8^5 

which successfully reduces the perturbative coefficients for every power of ui. Trying the 
same thing with the clover coefficient (145|) gives 

^piaq,MF^ !^ 1 1 + ^2^^ \Cf (-0.019865 + 0.079148CJ + 0.813321 u;2) 

+ iV, (0.015025 + 0.009617CU- 0.284786 cu^) |. (67) 

Again, mean field improvement works well. 

For the Symanzik action we calculate the integral in (163!) numerically, and get the 
result 

A;Sym = TT^ (0.732525- 11.394696 cu + 50.245225 cj^) . (68) 

The corresponding mean field improved expressions for k^ dSSD and csw (HH]) are 

^Sym,MF ^ [l + ^^ r<^ /QQOgQsa + 0.0500781 L^ - 0.471784 u;^)] , (69) 

8m5 '- ^ ^-' 

^Sym,MF ^ ^|l + gl^p \Cf (-0.0211635 + 0.115961 uj + 0.685247 w^) 

+ A/'c (0.013777 + 0.015905 cu- 0.321899 cu^) |. (70) 
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5.3 Choice of g\^p 

In this section we discuss the boosted couphng for S'f/(3), we have set Nc = 3, Cp = 4/3 
throughout. 

From higher order continuum calculations we know that g^-r^{fi) is a good expansion 
parameter if fi is close to the appropriate physical scale. On the other hand, series in the 
bare lattice coupling g'^{a) usually converge poorly. To understand this difference let us 
compare the two couplings. To one-loop order we have 

2bo (log -^ - log -^ ) = 26o log(a/i) + dg + Nfdf, (71) 



9^if^) 9^(0') V ^Ms «^ 



Mat 



where bo = (11— 2A^j/3)/(47r)^, and Nf is the number of flavors. The ratio of A parameters 
is thus given by 

Aiat _ fdg + Nfdf 



The coefficient dg is known for the plaquette and Symanzik gauge action |17j : 

rfPi^q = -0.4682 , df"^ = -0.2361 . (73) 

In Appendix C we show that dj is independent of the stout smearing parameter a;. There- 
fore, we can use the value for clover fermions computed in [18] 



For Nf = 3 this leads to 



rf/ = 0.0314917. (74) 



-^ =0.038 Plaquette, (75) 

-^ =0.289 Symanzik. (76) 

These ratios are far from 1, especially for the plaquette action, which explains the poor 
convergence of series in g'^{a). 

Now let us see what happens to the Lambda ratio if we make the popular choice of 
boosted coupling 

9mf = 71 ■ (77) 
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Upon inserting (164|) and (177j) in (ITTl) . we obtain 
which gives 



26o (log ^ - log -^ j = 2bo log(a/i) + d^ + iV^ rf^ + |^ , (7J 



For Nf = 3 the numerical values of this ratio are 



Ai^ _ ^^p / ci, + iV^d; + V37r^ ^ ^^g^ 



'/ 



-^ =0.702 Plaquette, (80) 

aMF 

-^s^ =2.459 Symanzik. (81) 

We see that mean field improvement drives Aiat towards Ajfg for both the plaquette and 
Symanzik gauge action, giving gl,jp ^ 9ws^ ^° ^^^^ 9mf appears to be a good expansion 
parameter in both cases. A perfect match is obtained for /i = 1/0. 702a (yU = 1 /2.459a) 
for the plaquette (Symanzik) action. 



6 Concluding remarks 

In the present paper we have computed the improvement coefficient csw and the additive 
mass renormalization/critical hopping parameter in one-loop perturbation theory for gen- 
eral stout parameter u performing a single smearing. To separate the effect of improving 
the gauge action from the effect of tuning the fermion action, we have done the calculation 
for both the plaquette action and the tree-level Symanzik gauge action. In addition we 
also present the 0{g^) corrections to the coefficients cngi and cd needed to 0{a) improve 
the quark fields in the most general case. 

We give mean field (tadpole) improved results for Kc and csw- For both the plaquette 
and the Symanzik action the boosted coupling gl^p turns out to be close to gj^, which 
makes gljp a good expansion parameter. We thus may expect that the perturbative series 
converges rapidly. 

For Nf = 3 fiavors of dynamical quarks it turns out that the one-loop improved 
Symanzik gauge action [2] coincides largely with its tree- level counterpart, with coefficients 
Co ~ 5/3, Ci ~ —1/12 and C2 ~ [T^. This makes the tree-level Symanzik action ([T]) 
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stand out against other improved gauge actions, at least from the perturbative point of 
view. SLiNC fermions represent a family of ultralocal, ultraviolet filtered clover fermions. 
While they share all prominent features of clover fermions, among them 0{a) improvement 
and flavor symmetry, they allow to further optimize the chiral properties of the action 
by tuning the fattening of the links. In our forthcoming simulations with Nf = 2 + 1 
and 2 + 1 + 1 flavors of dynamical quarks at realistic pion masses we shall employ this 
combination of gauge and fermion actions. 

Knowing the perturbative (asymptotic) value of csw: we can derive a closed expression 
for csw that covers the whole range of g^. We will do so in a subsequent paper employing 
the Schrodinger functional method. The one-loop coefficient Cgyy varies only slightly 
within the interval < a; < 0.2 for both the plaquette and Symanzik action. For u = 0.1, 
which is our favorite value, the tadpole improved one-loop coefficient becomes Cg{^ ~ 0, 
indicating that mean field approximation works well. The final result is cf{y ~ us/uq to 
a very good approximation for both gauge actions, where us is the average smeared link, 
found by measuring the smeared plaquette, and uq the average unsmeared link, found by 
measuring the unsmeared plaquette. 

This is to be compared with c;^ ~ ^h4) '^'^^^ fermions with no smearing. We therefore 
expect Csw to be a steeper function of g^ in the case of SLiNC fermions than for clover 
fermions. 

Stout link fattening reduces the additive mass renormalization considerably, with and 
without tadpole improvement, as expected. In fact, the critical hopping parameter k^. can 
be tuned to its continuum value of 1/8 for an appropriate choice of uj. We also confirm 
by early simulations with this action [20] that the spread of the negative eigenvalues 
is reduced by a factor of ~ 2 for u; = 0.1 and non-perturbative csw^ as compared to 
ordinary clover fermions. SLiNC fermions have many other appealing features as well. 
The renormalization factors of quark bilinear operators, for example, come out to be very 
close to unity, which hints at virtually continuum-like behavior. 
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Appendix A: Feynman rules 

In this Appendix we give the Feynman rules for quark-gluon vertices derived from action 
([3]) with single stout smeared gauge link variables in the Wilson part and general Wilson 
parameter r. The pieces in the vertices proportional to csw are denoted with V . They 
have been rederived using our notations and they agree with the Feynman rules given 
in [11]. In the vertices we denote the incoming/outgoing quark momenta by P\/p2- The 
incoming gluons are described by momenta /cj, Lorentz indices a,/?, 7 and color indices 
a,h,c = l,...,Nl-l. 

For the color matrices we have: 

N"^ - 1 
Cf = — , [^^ T''] = T'^T'' - T^T'' , {^^ T^} = T^T^ + T^T" (A.l) 

^,,,^^^,^^^,^^,^^^ ^,^,,^^^,^^^,^^,^^^ ^,^,,^^^,^^^,^^,^^_ 

We use the abbreviations 

s^(A;) = sin (^k^ , c^{k) = cos [^kf/j , s'^{k) = ^ sj(/c) , 

s\k,, k2) = J2 Sa«(^i + ^2) s^(fci - k2) = s\ki) - S\k2) . (A.2) 

For later use we give the bare massless quark propagator 

The structure of the Wilson quark-gluon vertices is 

W^i/.b2, Pi) = i c^{p2 + Pi) l^, + r Sf,{p2 + Pi) 

W2^{p2,Pi) = i Sf,{p2+pi) Jf,-r c^{p2+pi) . (A.4) 

Let us introduce the following functions to be useful in the definitions of the improved 
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vertices 

Va^,{k,Uj) = 6a^, + 4:UJVa^,{k) (A. 5) 

Vai^ik) = Sa{k) Sf,{k) -Sai,s'^{k) 

9appi{ki, /ca) = (5q/3 Ca{ki + /cs) s^(A;i - /ca) 

- ^a^ Ca(fc2) S/3(2/Ci + /Cg) + 5/3^ C/3(/i;i) Sc,(2/i;2 + /ci) (A. 6) 

Wa^i{kl,k2) = Sa{ki + k2) s^{ki-k2) -Sa^iS^{ki,k2) , (A. 7) 

Wai,{k,0) = Va^^{k) 

The qqg-vertex: V^{p2,Phki;csw,^) 

The qqg-vertex including stout smeared links and clover contribution is given by the 
expression (pi + ki = P2) 

V:{p2.Pi.kv,csw.^) = -^T'^^\4^(A;i,a7)iyi^(p2,Pi) + C5W'iC(A;i). (A.8) 

The stout smeared part shows the separation property mentioned in [6]. The clover part 
is given by 

i>(fci) = -2^T»^5^a«^c,(A;i)s^(2A;i) . (A.9) 

The qqgg-vertex: Vfp{p2,pi,ki,k2]Csw,'^) 
We define the qqgg-vertex as follows (pi + /ci + /ca = P2): 

K"^(P2,P1, A:i, A:^; cw, ^) = K!;'^ + KI^''^ + csw Kpih, ^2) • (A. 10) 

The stout smeared part is separated into two parts proportional to {T"-,T^} and [T°-,T^]. 
The anticommutator part shows the factorization property mentioned for two and four 
quark operators 

Vip'^ = lag^T-,T'}Y,Vc.,{k,,u)Vp^{k2,u)W2^{p2,Pi). (A.ll) 



23 



The commutator part is given by 



VIT = lag'[T'^,T'']ALuJ29c.fs,ih,h)W^,{p2,Pi). (A.12) 



Note that this part is proportional to u. The part oc csw has been used in the form 
V:'p{h,h) = i\g^[T\T'] \2a^p[2cp{h) c^{h) c^^h + k^) c^^h + h) (A.13) 
- Ca{ki) C/3(/c2)] +'5«/3 ^ o-^^ s„ ( fci + /C2 ) [S;,(2/i;2) " s^(2/ci) 



Both (1A.12I) and (1A.13I) vanish for tadpole diagrams along quark lines. 



The qqggg- vertex: V;x(i?2, Pi,ki,k2,h; csw, w) 



Pi^ 



We present that vertex contribution in the following form (pi + /ci + /i;2 + ^3 = P2) 

1 



V^g'i{p2,pi,ki,k2,k3;csw,^) = 7: a 9 x 



^/37 



6 



Yl \ ^i/^(P2,Pi) K%^,{ki, k2, ks) + cyclic perm. 



-QujW2^,{p2,Pi) T^a" Vaf,{ki) gp^^^{k2, ks) + cyclic perm. 
+ cswV:'^''^{k^,k2,k,). 



(A.14) 



Cyclic permutations have to be performed in the gluon momenta as well as in the color 
and Lorentz indices of the three gluons. Note that the general stout smeared part is 
proportional both to Wif^ and W2^- 

The coefficient F^^^^^{ki, /c2, /cs) is decomposed into its different color structures: 
F:%^{k,, k2, k,) = T,f /X(^i' ^2, ks) + T„t iCiki, k2, k,) - f%^{k,, k,, k2)) 



nbc \ p(3) 



+ [TT-^d'^"' a,,{kuk2,ks) 



N, 



(A.15) 
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where the /^.^^^^ are given as 



'a/37^ 



- Vaf^{kl, Uj) VfSi,{k2, Uj) V^i,{h, Uj) 



/SU(^i'^2,A;: 



^ VafMih, uj) VfSi,{k2, uj) 6^1, - - 5a^,5pi, V^f,{k3, uj) 



3) = 7^ yafiyKi,UJ) Vi3i_iyK2,UJ)Uyi_i— -Ua^iOi3ij,V^f^yKs,UJ) (A. 16) 

6uj6ap c^{ki - k2) cp{2k3 + ki + ^2) <^7/. + s^^ks) s^{k3 + 2ki) 6/3^ 

/a^V*^^l' ^25 ^3) = '2UJ 6f3^ (3 Wa^i{ki, /C2 + k^) + Va^i{ki + /C2 + /Cs)) (^q/S 

+ 12s^(A;i) Sa(/c2) s„(/c3) (s^(A;i + A;2 + /^s) f^a^ " ^»{ki + k2 + ks) Sfs^) 



.(3) 



The clover part of the qqggg- vertex is given by 



-abc 



V^iku ^2, ^3) = 27 i V^.Jki, k2, ks) + total perm. 



g I a/37 



(A.17) 



with 



~abc 

V ^p^{ki, k2, ks) 



T^'T'T'd^pd^-, 5Z ^"4 - ^ ^-(^1 + ^2 + k^) s^(2(A;i + ^2 + k^)) 
+ Caiki + k2 + ks) Cf,{ki + k2 + ks) Cf^^ks, - ki) s^{k2 

rparjibrjic i rjicrpbrjia 



(^al3 X 



(A.18) 



2 6f3^ Ca{ki + k2 + ks) cp{ki + /c2 + k^) c„(/c3 + k2) sp{ki) 

+ (5/3^ S/3(/C3 + k2) Ca{ki + 2/i;2) 

+ ^a-y Sa{ki + 2/^2 + k^) Cp{ki + k2 + ks) Cp{ks - ki) 



In (1A.17P the total permutation has to be performed in the gluon momenta, color and 
Lorentz indices. 

We only need this vertex for the gluon tadpole diagram of Fig. [H which simplifies the 
expressions. In the tadpole contribution to the vertex (1A.14P we denote the external gluon 
momentum by g = p2—pi, the color index of the gluon by a and the internal momenta by k 
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and —k. The color indices (6, c) of the remaining gluons forming the tadpole are summed 
up using the color diagonality 6^'^ of the gluon propagator, k is the gluon momentum in 
the tadpole loop. So the stout smeared tadpole contribution is defined from the general 
qqggg-vertex (explicitly symmetrized in the three gluons) as 



Af2-1 



6=1 



V:p^{p2,p„k) = Y.{KfM^PiA.k.-k) + cswV:'^^{v2,Vi.q.k,-k) 



\a''g^T'' J2 Wi^ip2,Pi)V^p^^{q,k) 



(A.19) 



N^-l 



CSW 



Yl KpM^Pi.q.k,-k). 



Using that definition we obtain for the stout smeared part 



Va^M^k) 



f{l) 



A^, 



p(2) 



f(2) 



(6 Cf - iV,) /^^,;^(g, k, -k) + -f [/^^^(A;, -k, q) - /^^^^(A;, g, -k) 
&,i-k, g, k) + f%,{-k, k, q)] + 4 Cf /X(g, k, -k) (A.20) 



+ (4 C^ - iVj ffl(k, -k, q) + /S,,(-A;, g, k) 



From that expression a convenient representation is found in the form 



Vap^^{q,k) = \ {6Cf- Nc) Vaf,iq,uj)Vf3f,{k,uj)V^f,{k,u;) 



+ 



2 

2u 



2 Sa^,Vp^{k, u) Vy^{k, u) - K^(g, uj) (Sp^ V^^{k, u) + S^^ Vf3^{k, u) ) 



3 (4 CiT - A^J C„;37m(?' ^) + Nc D^p^^{q, k) 



(A.21) 
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The structures Caf^^^ and Dap-y/j,, additionally contributing to O^u), are 

-ASy^Sfsip) ^a{k) [Saf3S^{p) S^{k) - Sa^,Sp{p) S^(A;) -Spi,Sa{p) Sa{k)] 

- 5a^5p^6^^ [2s\p) + 2s\k) -s\p + k)- s\p - k)] , 

(A.22) 
DafS'Yf.iq, k) = -3 5ay5p^ cp{p + k) Cy{p + k) - 3 5„^(5^^ c^(p - k) c^{p - k) 
+ 4 6pr^{5ap + 5/3^) Sq(p) s^(p) + 4 6a^i{Sp^ + 5^^) s^(A;) s^(/i;) 

- 2 (5,;.(5^^(5^^ [s^ip) + s2(A;)] + 6 5^^6p^ [2 c^ (p) c^(A;) - l] . 

Appendix B: Three-point function - universal part 



As discussed above, the universal part of the three-point function has the form ([3] 
when csw = 1 + (^{g"^). Therefore, it is sufficient to give only the one-loop result for 
Af|^(pi,P25 ?)• It is cast into the following form {q = P2 — Pi) 



Ki^{Pi,P2,q) = Fi{p^,p2)-f^ + F2{puP2) 



(B.l) 



+ [F3{Pl,P2) Pi,,, + F^{Pl,P2) P2.A f^ 
+ [F-,{PI,P2)P2,,, + ^6(Pl,P2)PlJ^2 • 

Due to the symmetries -^5(^1,^2) = -^3(^2,^1) and -^6(^1)^2) = -^4(^2,^1) we have four 
independent functions Fi{pi,p2) only. We represent them as follows: 



Fi(pi,p2) = 4CF^-^{12 + 2^-e)+2Q{C^S + N,p,.p2 + Cpq') 
CHl-a.f(4-0).og(|i 

+ V^ipuP2)l0g(^) +V,{p2,Pl)^Og(^ 



(B.2) 



F.(,„p.) = |(2iV.(6-0+c/^-^^^' 



A 



4A 



Pi.glog( ^j -p2.glogf -| 



,(B.3) 



i^3(Pl,P2) = C3f^ + 



P\ , 2Ar,e , e 



2A 

1 



q' 



8A 



+ - 



AN,i{p^.p2f+[2C^{QS + pl)- C^Pi.q)pl 

(B.4) 



V2ipuP2) log ( ^ j + V3{pi,p2) log ( -| 



27 



FiiPl,P2) 



C 



Pi.P2 2N,e e 



2A 



q' 



8A 



4(8Cp-iV,(4-0)(pi.P2)' 



[12CsS + 4Cepi + {C,+8CF{2 + 0)pt)Pi-P2+Crpipi (B.5) 






VA{pi,P2)\0g[-^\ + 1/5(^1,^2) log (^ 



The function Vj in front of the logarithms are found as follows 

VM.P2) = CH3 + 0-^(4-0+Ci^^, 



1^2(^1,^2) 



V3ipi,P2) 



4A 



(4C3 -C4 -4iV,0P2?' 

+ (12C3<S + 4iV,ePi-P2 + (C5 +8Ci.)g2)p2.g 
1 



4Ap^ 
^^4(^1,^2) = V2{p2,Pl 



[-AN,^pi.p2P2.qpl+ {-12C3Spi.q+ C^plq^) pj] 
1 



(B.6) 



4A 
1 



[-8 Cf{1 + Pi.q + (4 Cf{1 - 30 + iVc(5 - 00 P?] 



V^5(Pi,P2) = V2{puP2) + -r-r[i8Cp + N,i2-OOP2.q + {l + Oi'^CF + N,Opl]- 



4A 



We have introduced the kinematic functions 



A 

e 



{P1.P2) -P1P2, s 



pIpIq^ 



4A 



TT^ 



Sp 



P2-g+ VA 



pI 



Sp 



P2-g 



pI 



(B.7) 



1, 

+ 7^ lop 



P1-P2 



P1-P2 



A\ fq' 

A/ VP2 



with Sp(x) being the Spence function: 

Sp(x) = 



dy 



log(l - y) 
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The quantities Ci depend on the color factors and gauge parameter and have the values 

C2 = 8C^ + iV,(2 + (3-00), 

Cs = 4C^(l+0-iVc(4+(l-00, 

C, = 8C^(2 + 0-iVc(12 + (4-300, (B.8) 

C, = -iV,(4-(2 + 00, 



In order to express the one-loop result (IB. ID in terms of Spence functions, logarithms 
and rational functions of external momenta we have proceeded in two steps. First we have 
expanded all tensor integrals over the internal momentum into scalar three-point integrals 
times tensor functions of the external momenta [21] . Then we used recursion relations of 
Davydychev [22j to reduce these scalar three-point integrals into scalar two-point integrals 
and G. 



Appendix C: ct;- independence of df 

We find df, the coefficient which tells us the fermionic shift in Aiat, by calculating the 
massless quark vacuum polarization in a gluon with a^g^ -C 1: 



^afsiQ-.Csw,^) 



- Nf 



Nf 



(2^ 
d^k 

(M 



Tr \V^{q + k, k, q; Csw, u;)S{k)V^{k, q + k, -q; Csw, u;)S{k + q)] 



rTr 



r{a,b} 



V/3' (^' ^' ^' ~^' ^sw, u;)S{k) 



(C.l) 



The quark propagator S and the vertices V are defined in Appendix A, the trace here is 
over both spin and color. The corresponding one-loop diagrams are shown in Fig. HI 

In the required limit of small a^q^ we can expand in g^ and drop any terms 0{a^q'^). 
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(a) 



(b) 



Figure 4: One-loop quark vacuum polarization diagrams. 



We then get 



Uf^{q;csw,u;) = W^^^iq^Csw^O) - 2lj Nj S'^'g^a^ x 



J2 i^o^q, - q^5^,) j ^Tr [W,,{k, k)S{k)W,p{k, k)S{k)] 



+ a {qaQp - q^Sap) 



Yl (^/39m - q^h^^) 



+ a [q^^qp - q^5ap) 



d^k 

(2^ 



Ti[W2p{k,k)S{k)] 



(C.2) 



d^k 

(2^ 



i:i[w^^{k,k)S{k)w^,{k,k)S{k)] 



d^k 

(2^ 



2„4n 



Ti[W2^{k,k)S{k)] )+0{a'q 



where n^^(g; Csw^ 0) is the vacuum polarization tensor with no smearing, Wi and W2 are 
the Wilson quark gluon vertices defined in (lA.4p . and the trace is now only over the spin 
index. All uj'^ terms have dropped out because they first appear at 0{a?q'^). Calculating 
n^^(g; csw^ 0) in one loop for csw = 1 leads to the value of df given in Eq. fITil) . 



From power counting we would at first expect the integrals oc ci; in (]C.2p to have values 
proportional to 1/a^ or 1/a'^, and to make a finite contribution to df. However we show 
now that there is a perfect cancellation between the continuum-like diagram Fig. Hl^a) 
(the integrals involving Wi) and the tadpole contribution Fig. |l](b) (those with W2)- To 
do this we use the identities 

^ -S{k) = -S{k)Wi^{k,k)S{k), (C.3) 



_d_ 

dk„ 



dkf, 
Wi,{k,k) 



a6^,uW2^{k,k) 



(C.4) 
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which follow immediately from the definitions. Eq. (1C.2I) becomes 



{ E {qaq, - q'S^,) J ^ ^Tr [W,,{k, k)S{k)] (C.5) 

+ E iw, - q%,) I (01 ^Tr [W,,{k, k)S{k)] I + 0{aY) • 



The integrals are now zero because Wi and S are periodic, 

fvr/a o 

-7r/a 



r " rfA;„T|-Tr [W,^{k, k)S{k)] = Tr [l^i^(fc, k)S{k)] '" ''^'^ = . (C.6) 

J-Tr/a C>'^a ka=-'iT/a 



Thus we have proved that the vacuum polarization is independent of smearing the one-link 
part of the fermion action, 

Ii%{q- csw, uj) = n^^^(g; csw, 0) + 0{aY) (C.7) 

which implies that df depends on r and csw, but not on u. 
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